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Abstract 

.•^ . A general algebraic method of quantum corrections evaluation is presented. Quantum 

l/^ \ corrections to a few classical solutions (kinks and periodic) of Ginzburg-Landau (phi- 

CN ■ in-quadro) and Sin-Gordon models are calculated in arbitrary dimensions. The Green 
function for heat equation with a soliton potential is constructed by means of Laplace 

^T I transformation theory and Hermit equation for the Green function transform. The gen- 

QQ ' eralized zeta-function is used to evaluate the functional integral and quantum corrections 

^D ■ to mass in quasiclassical approximation. 

>: 

^j ' 1 Introduction. 

P.. 

1.1 General remarks. 

We consider one-dimensional field theory, based on nonlinear Klein-Gordon equations, arising, 
for example of Sine-Gordon (SG) case, in kink models for crystal structure dislocations [T]. In 
[2] the authors study the diffusion of kinks. The Ginzburg-Landau (GL) model is very popular 
in different aspects of solid state physics, e.g. for magnetics [3]. 

1.2 Feynmann quantization of a classical field. 

An attention to Feynmann quantization formalism of a classical field was recently attracted in 
connection with a link to a SUSY quasiclasic quantization condition [1] suggested in [5]; see, 
however, [H]. 



A class of nonlinear Klein-Gordon equations in the case of static one-dimensional solutions 
is reduced to 

0"-y'(0) =O,0 = 0(x),a;Gi?. (1.1) 

Suppose the potential V{(f)) is twice continuously differentiable; it guarantees existence and 
uniqueness of the equations correspondent to (11. ip Cauchy problem solution. The first integral 
of (11.11) is given by 

W = ^{^''-V{^), (1.2) 

where W is the integration constant. The equation (II. 2p is ordinary first-order differential 
equation with separated variables. As the phase method shows, solutions of this equations 
belong to the following families: constant, periodic, separatrix and the so-called "passing" one 
[7|. In the case of phi-in-quadro (GL, non-integrable) model, the "potential" is 

while in the case of Sin-Gordon (SG) model it is 

, , 2m'^ , , 1 /Sq , , , , 

K,(0) = ^(1 + cos(-y y^))- (1-3) 

We modified this last model to fit it with the first one at small values of the constant g, namely 

^.. = y^ + ^.^ + 0(/). 

In the papers of V.Konoplich [8] quantum corrections to a few classical solutions by means 
of Riemann zeta-function are calculated in dimensions d> \. Most interesting of them are the 
corrections to the kink - the separatrix solution of the field 0^ (GL) model. The method of [8] is 
rather complicated and it is desired to simplify it, that was the main target of our previous note 
[9]. We applied the Darboux transformations technique with some nontrivial details missed in 
[8]. The suggested approach open new possibilities; for example it allows to show the way 
to calculate the quantum corrections to matrix models of similar structure , Q-balls [10] and 
periodic solutions of the models. The last problem is posed in the review [TJ. 

The approximate quantum corrections to the solutions of the equation (ll.ip are obtained 
via the Feynmann functional integral method evaluated by the stationary phase analog [H]. It 
gives the following relation 

exp[-^] ^ -^, (1.4) 



where Squ denotes quantum action, corresponding the potential ^(0), A - some quantity de- 
termined by the vacuum state at ^(0) = 0, and detD is the determinant of the operator 

D = -dl-^,-.V'\<\>{x)). (1.5) 

The variable y G R^^^ stands for the transverse variables on which the solution (^{x) does not 
depend. The operator D appears while the second variational derivative of the quantum action 
functional (which enter the Feynmann trajectory integral) is evaluated. For the vacuum action 



Syac the relation of the form (11.41) is vahd if Squ is changed to Syac and D is replaced by the 
"vacuum state" operator Dq = —d^ — A.^. Then, the quantum correction 

is obtained by the mentioned twice use of the formula (11. 4p as 

. ^ h^ . detD . , „. 

Hence, the problem of determination of the quantum correction is reduced to one of evaluation 
of the determinants D and Dq ratio. The methodic of the evaluation will be presented in the 
following section. 

1.3 The generalized Riemann zeta-function and Green function of 
heat equation. 

The generalized zeta-function appears in many problems of quantum mechanics and quantum 
field theories which use the Lagrangian £ = (50)^/2 — V{(j)) and it is necessary to calculate a 
Feynmann functional integral in the quasiclassical approximation. 

The scheme is following. Let {A„} = S* be a set of all eigenvalues of a linear operator L, 
then, logarithm of the operator determinant is represented by the formal sum over this set 

In(detL) = ^lnA„, (1.8) 

A„es 

where the sum in the r.h.s. is formal one, S means the operator's spectrum. The generalized 
Riemann zeta-function Cl{s) is defined by the equality 

Cl{s) = Y. X-'. (1.9) 

A„e5 

This definition should be interpreted as analytic continuation to the complex plane of s from 
the half plane Res > a in which the sum in (11.91) converges. Differentiating the relation (11.91) 
with respect to s at the point s = yields 

ln(detL)=Ci(0). (1.10) 

The generalized zeta-function (I4.7ip admits the representation via the diagonal qd of a 
Green function of the operator dt + L. Such representation is obtained as follows. 

Let r G M be the set of independent variables of the operator L; particularly, the operator 
D of (11.51) depends on r = (x, y) G M x M'^"^ then, at t > 0, 

idt + L)g{t,T,To)=5{r-Vo) (1.11) 

and 

^(t,r,ro) = 0, t<0. 

There is a representation in terms of the formal sun 

9D{t, r, ro) = Y^ exp[-Ant]^n(r)C(ro), (1.12) 



where the normahzed eigenfunctions ^ra(r) correspond to eigenvalues A„ of the operator L. Put 
L = D. Let us introduce the function 

Init) = / drgnit,r,r) = ^exp[-Ant], (1.13) 

■^ n 

that follows from (11.121) and normalization. The integration in fll.l3p is performed either 
along the axis in the kink case or over the period in a case of periodic solutions. The Mellin 
transformation of fll.lSp yields in 



hoo 



r(s) Jo 

where the T{s) is the Euler Gamma-function. 

The generalized zeta- function, defined by the relations (I1.13fl.l4p . will be referred as the 
zeta-function of the operator D. 

From the relation (11.131) for the function 7D(t) it follows an important property of multi- 
plicity: if the operator D is a sum of two differential operators D = Di + D2, which depend on 
different variables, the following equality holds 

lD{t)=lDAt)lDM- (1-15) 

We will need the value of the function 7D(t) for the vacuum state, when the operator 
D = Dq = —A is equal to the d-dimensional Laplacian. In this case the formal sum in the r.h.s 
of (11.131) goes to d-dimensional Poisson integral 

^^oit) = tA^ / dkexp{-\k\H) = {Antr'^/\ (1.16) 



More generally, for a constant potential z/. 



l-d i (S H TT' I d-l~2s 



7D.(t) = iAn)^'-^^^^^u^^^. (1.17) 

The basic relation (I1.7P points to a necessity of evaluation of the determinants of the oper- 
ators 

D = Do + u{x), Do = -dl-Ay + \ (1.18) 

where A is a positive number and x G -R is one of variables, while y G Rd-i is a set of other 
variables. The operator Ay is the Laplace operator in d-1 dimensions, u(x) is one-dimensional 
potential that is defined by the condition 

V^"(0o(a;)) = A + M(x), (1.19) 

where (poix) is the classical static solution of the equation of motion. 

A quantum correction to the action in one-loop approximation for the classical solution 0(x) 
is calculated via zeta-function by the formula 

Ae = -C^(0)/2, (1.20) 



where 

r-OO 



/■oo 

CDis)=M'-' / ^it)f-'dt/Tis); (1.21) 



here r(s) is the Euler gamma function and M is a mass scale. 

The function 7(t) in the MeUin integral fll.2ip is expressed via the Green functions difference 
G{x, y, t; Xq, yo, to) and Go{x, y, t; Xq, yo, to) oi dt + D and dt + Do in the following way; let 

g{x, t) = G{x, y, t; xq, yo, 0) - Goix, y, t; xq, yo, 0), 

due to the translational invariance along y of operators D and Do the function g does not 
depend on y; the contraction of Go is necessary for deleting of ultraviolet divergence. 

2 Static solutions of the SG and (p^ model 

Starting with SG model we integrate the equation (11. ip with the potential (11.30 arriving at the 
first-order differential equation with the parameter W (II. 2p . 

yf = ^^(1 + cos{-\Uv)) + 2W. (2.22) 

Sg m \ 2 

The solutions are restricted and hence have the direct physical relevance, if 

4m^ 



3^ 



<W <0, (2.23) 



that follows from phase plane analysis. li W = the nontrivial solutions are interpreted as 
kink and antikink 

[Y [Y 

(p(x) = ±\ — arcsintanh(ma;) (mo(i$), $ = 2rmT\ — , (2.24) 

V 35- V 35- 

while at the interval 

<W <0, 

the solution of (12.221) yields a periodic function expressed via elliptic Jacobi function. To find 
it one plug ip = ±2mj- arcsinz, then the equation (I2.22p goes to 

(.')^=m^(l-.^)(l + |^). (2.25) 

The solution of (12.251) at the interval (I2.23P is given by 

z = ksn{mx\ k), (2.26) 

where 

is the module of the elliptic function. Hence 

A(k^ - l)m^ 
W = - — — ^ — . 2.28 

3^ 

5 



Finally 

2 

if = ±2m — arcsin ksn(mx;k) (mod^). (2.29) 

The class of restricted solutions contains also 

i777 

^ = 0(mod<l>);W = (2.30) 

3^ 

and 

^ = ±nmJ — (mod<l>);W = 0. (2.31) 

V 3^ 

Other static solutions are obtained by shifts 

that follows from Klein-Gordon equation invariance and SG equation potential periodicity. 

In the case of static solutions of the ip^ model the potential and the Lagrangian are deter- 
mined by the formulas 

V{v) = gif^A - mVV2, 

therefore the equation of motion has the form 

<p"{x) + mV - 9V^ = 0, (2.32) 

that yields (11.21) in the form 

{^r = k^'-—? + 2W. (2.33) 

^ 9 

Its restricted solutions, as it follows from phase plane analysis, exist if 

2 
771 

<W <0. (2.34) 

The separatrix (W=0) solution of (12.331) is the kink/antikink 

/2 m 

-6tanh(6x)), b = ^. (2.35) 

9 V2 

While inside the interval (12.341) the equation (I2.33p is expressed in terms of the elliptic 

Jacobi sinus 

[2 rn 

ip = ±J-kbsn{bx;k), b= < k < 1. (2.36) 

The constant W is given by 

The family of the restricted solutions contains also the constant vacuum ones (W=0) 

^ = ±^. (2.38) 

6 



After the substitution of (12.351) into (11.191) we obtain the following potential u(x) : 

u{x) = -6b'^/ch^{bx), (2.39) 

with the meaning of the constant b = m/^/2. As a result the two- level reflectionless potential 
of one-dimensional Schrodinger equation —d"^ + u{x) appears. Eigenvalues and the normalized 
eigenfunctions of which are correspondingly (its numeration is chosen from above to lowercase). 



Ai = —6^, ^i{x) = a/36/2 smh{bx) / cosh"^ (bx); 
A2 = -46^, ip^[x) = V3b/2cosh{bx). 



3 The energy of classic static solutions of SG and 0^ 
models 

Let us evaluate the energy of the nontrivial static solutions of both models via the energy 
density definitions 

POO 

eix)= / ii<t>T/^ + Vi^))dx, (3.40) 



for kinks and, in a case of periodic solutions, 



E = 2 I e{x)dx, (3.41) 

Jo 



the constant "1" is the period of the solution. 
For SG kink: 

E, = ^, (3.42) 

9 



^m^,_ ,2\ 



for a periodic soliton 

where K(k),E(k) - complete elliptic Legendre integrals. 



Ep = [{l-e)K + 2E], (3.43) 



4 The generalized zeta-function as analytic function. 

Let us consider the Laplace transform of the one-dimensional Green function defined in the 
Seel by the relations fll.liyi.12p 

gL{t, X, xo) = -— hip, X, xo)eP\ (4.44) 

2^ Ji 

where the integration is performed along the contour in complex p-plane that contain all the 
singularities of the integrand . The integrand function gL{p,x,xo) may be considered as the 
Green function of the spectral equation 

{L + p)gL{p,x,xo) = 6{x-xo). (4.45) 



The corresponding transform of 7(t) is denoted as 

l{p)= 9L{p,x,x)dx, (4.46) 

the inverse transform of (I4.46P looks as (14.441) with the same contour of integration 

7L(t) = f^{p)eP'dp. (4.47) 



The inverse of the formula ( 14.46P may be considered as a base for analytic continuations. It 
is known from the general properties Laplace transform of a distribution that are zero at t < 
[H], so that the functions (14.461) and gL{p,x,xo) are analytic function at Rep > ctq- 

It also shows that no necessity to solve the equation for the Green function (11.111) . it is 
enough to solve the spectral problem (I4.45P 

In a spirit of Hermit approach, see, e.g. [15] , the function giip, x, x) is a solution of bilinear 
equation 

2GG" - (G'Y - A{u{x) - p)G^ + 1 = 0, G(p,x) = gL{p,x,x) (4.48) 

which in a case of reflectionless and finite-gap solutions is solved more effectively than (14.451) . 
In the case of the periodic solution of the phi-in-quadro model the potential has a " cnoidal" 
form 

u{x) = -Qk%^cn^{hx] k) + (hk'^ - l)h^. (4.49) 

This potential differs from second Lames equation potential by the constant hence its spectrum 
is two-gap one. 

It is possible to cover all necessary classes of solutions of both models (cases A,B for 0^ and 
C,D for SG) via the universal representation by means of polynomials (in p) P, Q 



where 

for kinks A,C, and 



G{p, x) = P{p, z)/2^Q{p,z), (4.50) 

z = sech'^ipx) 

z = cn^ipx] k) 

for the periodic B,D. 

In the 0^ (A,B case)the function P(p) is linear in p polynomial while Q(p) is one of the 
third order. Plugging (14.501) yields 

h\p{z){2PP" - {P'f) + p\z)PP') -{p + u{z))P^ + Q = 0, (4.51) 

the primes denote derivatives with respect to z, while 

p = { z'^{l- z), cases A,B; (4.52) 

z{l- z){l-P + k'^z), cases G,D. 

b'^{l-2z),case (A) 
,, , b'^(2k'^ -l-2k'^z),case (B) ,, ^^, 

^ 2b [2 — 3z), case (6) 
b'^{5k'^-6k'^z),case {D) 

8 



Let us start with the case (A). The form of the polynomial is determined from well-known 
facts of the reflectionless potentials theory. 

Q=p^{p + b^). (4.54) 

The case (B) 

P = p + Pi{z), Q = p^ + q2p'^ + qip + qom (4.55) 

includes (A) in a sense that for the case q2 = b'^, qi = 0, go = 0. 
Substituting (14.551) into (14.511) gives for each power of p = 0, 1, 2 

-2Pi(z)-n(z) + g2 = 0, 
b\2p{z)Pi' + p\z)P[) - Pf - 2u{z)Pi + gi. (4.56) 

b\p{z){2PiP[' - P2) + p'{z)P,P[ - u{z)Pl + go = 0. 

respectively. In the case (A) it gives 

Pi = b'^z. (4.57) 

For the case (B) we begin from 

Pi = k'^b'^z, (4.58) 

that is the result of substitution of Pi from the first equation of ( 14.56P into the second one. 
Next, from the third equation yields 

Q = p^p'j^(l-k^)b''){p-k%''), (4.59) 

with the simple roots pi. 

Going to the cases (C,D), generally 

P = p'^ + Pi{z)p + P2{z), Q=p^ + q4P^ + q3P^ + q2P^ + qiP + qo, (4.60) 

while in the particular case of (C) 

go = 0, gi = 0, ga = 366^ gs = 336^ gs = lOfel 

A substitution of (I4.60p into the Hermit equation splits in the system 

-2Pi - u + g2 = 0, 

-2P2 - Pf - 2uP, + 62(2pPi" + p'P{) + gg, 

b\pi2P^ + 2PiPf - (Pif) + p'iP^ + P,P[)) - 2P1P2) - «(2P2 + P2) + g2 = 0, (4.61) 

b\2p{2P['P2 + P[P^ + PiP^O + p'(PiP^) + Pi'P2)) - Pi - 2uP,P2 + gi, 

b'{p{2P2Pi'-Pi'))+p'Pi + qo. 

The arguments in (14.611) are omitted. Solving the system yields for the case D 

Pi{z) =az + (3, P2{z) = az"^ + bz + c, (4.62) 

where 



a = 3k^b\ 
2 



^ (l-5fc2)b^-(?4 



a = 9k%\ (4.63) 



2 ' 

4g3-9|+2(l-5fc^)fe^g4+3(l-6fc^+21A:-^)fe-^ 

8 



are functions only of k, b, where 



g2 = -9(1 - 3P - 3A;4 - A;6)66^ q^ = 3(1 + p + k^)b\ q^ = 5(1 + e)b^. ^^'^^^ 



Finally, 



i=5 



Q = l[{p-Pi), (4.65) 



i=l 



where the polynomial Q the simple roots Pi are ordered so that (Ojkjl) 



- (2VI - P + A;4 + 1 + ^2)^2 ^ _3^2 ^ _3^2^2 ^ Q ^ -(2VI - P + A;4 - 1 - e)b'^. (4.66) 
Let us pick up the expressions determining 7(p): 

So we need three integrals over the period K. 

/g dx, /q 2;(ix, /g 2;2(ix. (4.68) 

Let us go to the variable z, dz = d{cn'^{x)) = —2cn{x)cn{x)dn{x)dx = 2sn{x), a bit more 
convenient to put y = 1 — z,dy = —dz 

Jo ^^ = Jo y{l-y){l-k2y) = 27^(^)5 
Jo ■^''^ v^' ^)'^^ ^ 2 Jo y{l-y){l-k2y) ^ P ' (4.69) 

/g^ .^^(x; fc)rfy = /g^ ,(,.;;(t,.,) = 3f((2 + A;2)i^(A:) - 2(1 + P)E{k)). 



hence 

7(p) = 2K{k)^ + a(-^(^^^ + (/? - a)2i^(fc)^+ 
^^ i2+k^)Kik)-2ii+k^)Eik) ^ ^^ _ 2a)^(^^ + (a - 6 + c)2ir(/t))^. 

or, finally 



(4.70) 



Cis) = M'^ J {j ^ip)eP'dp)f~'dt/Tis); (4.71) 

which is expressed via integrals 

L '$dP (4.72) 

by a contour that contains all branch points of the integrands. 

5 SG kinks 

The results of the previous sections allow to evaluate corrections to actions for all four (A,B,C,D) 
cases. The results for cjp' model kinks are well-known, see, e.g. [9j, where a table for the 
dimensions d=l,2,3,4 is listed. These results fit the case B, which hence was strictly verified. 

10 



Let us present the formulas for the kinks of the SG niodeL The substitution of the expres- 
sions from fl4.50|4.65|4.55|4.57l) yields a divergence of the integral for (14.671) 7(p). To regularize 
the Green function let us divide itd Laplace transform into two parts as 

G{p,x) = G,{p) + Gk{p,x), (5.73) 

where the part Gc{p) is the Green function diagonal (the solution of (14.481) for a constant 
potential: 

Ge = — j^ . (5.74) 

The kink part is easily constructed via (14. 50f4.65f4. 55114. 571) : 

G, = ^!£^£i!M. (5.75) 

1p\/p + 6^ 

The representation (I5.73P results in two contributions of the quantum corrections to the kink 
(antikink) mass. 

The first one coinsides with (15.761) for v = m^ . 

7^oW = (4vr)V£if±M^^---, (5.76) 

r(s) 

and the second one is obtained by the general scheme, namely 

/oo 
Gk{p,x)dx. (5.77) 

-oo 

Plugging (I5.75P into (I5.77p yields 

7.(p) = /— - ■ (5.78) 

PyP + ^ 

The value of the integral J_ sech'^{bx) = 2/6 is taken into account, the corresponding function 
(ILlSp or, more exactly (14.470 will be denoted j^, the index R label the result of a renormal- 
ization provided by the division (I5.73p . The function may be found directly from a table ([17]) 
but we would explain the result as an example for further development of the renormalization 
procedure. Note that the cut for the radical ^/p + IP is made along the Re p -axis from — cx) to 
—6^, the branch with a/p + 6^ = i^/\p + lp\ on the upper and \/p + b^ = —i\/\p + 6^| on the 
lower bounds of the cut is chosen. After the deformation of the integral contour one has for 

The formula (15.791) gives an expansion of the resolvent of the operator dt + Lhj the operator L 
spectrum. The direct substitution of (I5.79P into (11.141) leads to a divergent integrals. However 
in this case a renormalization is not necessary, the integral in the r.h.s. of (I5.79P is expressed 
in terms of the error function 

Erfiz) = ^ r e-^'dr. 
11 




Figure 1: d=l, dependence of the correction Ae to mass on the parameter m 



namely, after some change of variables 



7fc(t) = ErfibVt) = ^ f expl-bHr^T. 

V^r Jo 



(5.80) 



The same result gives [T7J. The Mellin transform of this representation gives the following 
expression of zeta function of the operator L. 



C."(«) 



26-2^ r(. + 1/2) f _,^_, 



Tis] 



-2s 



'dr 



r(s + 1/2) 



(5.81) 



,0 0r r(s + l) 

The integral in ( ]5.8ip converged only at Res < but gives the analytical continuation for all 
s E C excluding the poles of r(s + 1). 

Substituting the result (15.811) into one arising from (I1.15P with account of (I5.76P yields 

r(s + 1 - d/2) 



CE{s) = -4(47r)2m' 



d-l-2s 



{2s + l-d)T{s)' 



(5.82) 



The condition 6 = m is taken into account. Differentiation of (I5.82p by s at the point s = 
give the desired correction (11.201) 



-2(4:71)2 171' 



d-2s-i n'-¥+v 



1 <C£is)) _ 

2 ds 

ris)i2s-d-,ir ((^ -2^-1) (Psi Hd + s + 1) - 2 Inm - Psi (s)) + 2) 

(5.83) 

next we plot the dependence of the correction -^^-^^ — - on m (Fig [3]). 



We would remind about the choice of the constant g=2 in the last sections. 
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Figure 2: d=2, dependence of the correction Ae to mass on the parameter m 




Figure 3: d=3, dependence of the correction Ae to mass on the parameter m 
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6 Conclusion 

In the next paper we admit that the potential u(x) from f ll.lQp has the form of n-level reflec- 
tionless potential 

u{x) = -n{n + I)}?/ cosh^(te). (6.84) 

with eigenvalues A^ = —m'^b'^, m = 1, ...,n. We would note for further generalization means 
that this function may be considered as degenerate limit of n-gap Lame potential of Hill equa- 
tion. Such potentials correspond to higher solitonic models. 

The kink case corresponds n = 2, X = n'^b'^. The quantum correction to its action will be 
calculated in sec. 4 from general formulas of sec. 3. 

Let us note that all results related to the scalar nonlinear Klein-Gordon equation may be 
applied directly to many-component model 0(x) G §^ with account of SO(r?7,) symmetry. The 
equation fll.ip takes the form 

-r+v'm\)Aj = o- (6.85) 



The scalar operator D is defined by fll.Sp goes to the matrix one 

D = i-dl -Ay + ^^'(11011)^]/™ + [(^^"(11011)) - nil0ll)^]- (6.86) 



The technique developed in this paper is transported to quantum corrections for periodic static 
solutions of 0^ model: 

, , km [2 , mx , x , , ^ 

M^) = Y^,^j-sn{^===,kl 0<k<l; (6.87) 

where A; is a module of the elliptic function. When k = 1 the formula (16.870 goes to one for kink 
(I2.35p . The substitution of (16.870 into (11.190 yields in two-gap Lame potential that is embedded 
in Darboux Transformations theory by chain representation [TU] that give a possibility to derive 
the Green function analogue for this case. The results will be published elsewhere. Some recent 
papers open new field for applications [201 [21] . 
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